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Abstract. We quantify the degree of lopsidedness for a sample of 149 galaxies observed in the near-infrared from 
the OSUBGS sample, and try to explain the physical origin of the observed disk lopsidedness. We confirm previous 
studies, but for a larger sample, that a large fraction of galaxies have significant lopsidedness in their stellar disks, 
measured as the Fourier amplitude of the m = 1 component normalised to the average or m = 0 component in 
the surface density. Late-type galaxies are found to be more lopsided, while the presence of m = 2 spiral arms 
and bars is correlated with disk lopsidedness. We also show that the m = 1 amplitude is uncorrelated with the 
presence of companions. Numerical simulations were carried out to study the generation of m = 1 via different 
processes: galaxy tidal encounters, galaxy mergers, and external gas accretion with subsequent star formation. 
These simulations show that galaxy interactions and mergers can trigger strong lopsidedness, but do not explain 
several independent statistical properties of observed galaxies. To explain all the observational results, it is required 
that a large fraction of lopsidedness results from cosmological accretion of gas on galactic disks, which can create 
strongly lopsided disks when this accretion is asymmetrical enough. 
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1. Introduction 

It has been known for a long time that the gaseous com¬ 
ponent (HI) in galaxies is often strongly asymmetric and 
lopsided. Baldwin et al. (1980) have first studied this re¬ 
markable phenomenon in detail and proposed that the 
perturbations are m = 1 waves built from off-centered 
elliptical orbits that may persist a long time against dif¬ 
ferential precession. However, they have derived a lifetime 
that is still not sufficient to explain the high frequency of 
lopsidedness in neutral gas disks. Richter & Sancisi (1994) 
have compiled a sample of 1700 galaxies from the litera¬ 
ture and noticed on resolved maps of nearby galaxies (like 
M101) that there is a correspondence between lopsided¬ 
ness and the global HI velocity profile of the galaxy. It 
is therefore possible to study the asymmetries directly on 
the global HI spectrum, for a much larger sample. They 
derive a lower limit of 50% for the fraction of lopsided 
galaxies, and conclude that HI asymmetries in disk galax¬ 
ies are the rule rather than the exception. Lopsidedness is 
more frequent in late-type galaxies, where Matthews et al. 
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(1998) found a frequency of 77% of Hi-distorted profiles. 
However, it is difficult to get a precise quantitative indica¬ 
tor of the degree of lopsidedness with HI profiles without 
any spatial resolution. 

The ubiquity of lopsidedness may not be the privilege 
of HI disks. The stellar disk potential, traced by near- 
infrared light, can be strongly lopsided, even in isolated 
galaxies, like the spectacular case of NGC 1637 (Block 
et al. 1994). The stellar light is definitely affected by the 
phenomenon, as shown by the studies of Rix & Zaritsky 
(1995) and Zaritsky & Rix (1997) with a total sample of 
60 galaxies. They computed the mn = 1 Fourier ampli¬ 
tude of the density, and found that at least one third of 
these galaxies are significantly lopsided. Kornreich et al. 
(1998) quantify the asymmetry by comparing the optical 
light contained within trapezoidal sectors of galactic disks, 
and find also that 30% of galaxies are very lopsided. In 
addition to the spatial asymmetry discussed above, some 
galaxies also display kinematical lopsidedness though this 
has been less studied (Schoenmakers et al. 1997, Swaters 
et al. 1999). 

To explain the high frequency of lopsidedness, two 
kinds of arguments have been proposed. First the per- 
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turbation might be longer-lived than previously thought, 
either because the winding out by differential precession 
is quite long in the outer parts of galaxies (Baldwin et al. 
1980), because of weakly damped modes (Weinberg 1994), 
because of the amplification of density waves (Shu et al. 
1990, Syer & Tremaine 1996), or alternatively that the 
disk is distorted in a lopsided halo (Jog 1997, 2002) or 
else off-centered with respect to a massive halo (Levine & 
Sparke 1998, Noordermeer et al. 2001). Second, the per¬ 
turbations could be forced by tidal interactions (Kornreich 
et al. 2002) or minor mergers (Zaritsky & Rix 1997). In 
the last cases, the perturbations could be short-lived, and 
either lopsided galaxies should be observed with compan¬ 
ions or the merger frequency is very high. For instance, 
Walker et al. (1996) have shown through N-body simula¬ 
tions that a minor merger can result in a lopsided disk with 
perturbations lasting at least 1 Gyr. But until now, there 
has been no observed strong correlation between the pres¬ 
ence of companions and lopsidedness (Wilcots & Prescott 
2004). Thus, what generates the large fraction of lopsid¬ 
edness in the disks of spiral galaxies remains unknown. 

In this paper, we measure the lopsidedness for a larger 
statistical sample of galaxies, and try to understand the 
physical origin of the observed disk lopsidedness. In Sect.|3 
we quantify the m = 1 Fourier amplitude of the surface 
density in a sample of 149 galaxies with available near- 
infrared images from the Ohio State University Bright 
Galaxy Survey (OSUBGS, Eskridge et al. 2002). Since the 
galaxies chosen are not lopsided a priori, this sample is 
expected to give a definitive measure of lopsidedness in 
typical disk galaxies. We also measure the HI asymme¬ 
try for the present sample of galaxies. The dependence of 
the degree of lopsidedness versus other physical parame¬ 
ters (presence of companions, bars and spiral arms, and 
Hubble type) is studied in Sect. 3. In order to understand 
the origin of the disk lopsidedness, we examine various sce¬ 
narios via numerical simulations: distant tidal encounters, 
minor mergers, and gas accretion (Sect. 4). We show that 
galaxy interactions and mergers are not enough to explain 
all the properties of the observed galaxies. Cosmological 
accretion of gas seems to be responsible for a large frac¬ 
tion of lopsided disks. Observations of lopsided disks im¬ 
ply that galaxies have accreted large amounts of gas in 
past Gyrs, as already expected from the study of barred 
galaxies (Bournaud & Combes 2002, Block et al. 2002). 
These results are discussed and we give our conclusions in 
Sect. 5. 

2. Analysis of observations 

2.1. The OSUBGS sample 

To quantify the amplitude of the m = 1 mode both 
in the disk surface density and in the potential, we 
use near-infrared images from the Ohio State University 
Bright Galaxy Survey (OSUBGS, Eskridge et al. 2002). 
This is a magnitude-limited sample drawn from the 
Third Reference Catalog of Bright Galaxies (RC3, see 


Vaucouleurs et al. 1991). For a spiral galaxy to be included 
in OSUBGS sample, its B magnitude has to be brighter 
than 12 m , and its diameter less than 6'. According to the 
RC3, the catalog is at least 95% complete for galaxies this 
bright. We use the same techniques for deriving the vari¬ 
ous m-Fourier amplitudes as in Block et al. (2002). 

From the 174 galaxies in the sample, 171 are appropri¬ 
ate for studies of surface density distribution. We include 
only galaxies that are inclined up to 70° with respect to 
the line of sight, which gives a total of 149 galaxies. The 
choice of this inclination cut-off is justified in Sect. 2.3. 

2.2. The bulge/disk decomposition and deprojection 

The bulge and disk components are decoupled through 
radial profile Fitting; the bulge in the center is fitted by an 
r 1 / 4 de Vaucouleurs profile, and the disk by an exponential 
profile. 

To deproject the galaxies, a particular 3D shape has 
to be assumed for each component. We assume the disk is 
flat, with constant scale-height (which will be varied next). 
Bulge components are likely to be flattened, but less than 
disks, and the flattening of the bulge is quite difficult to 
derive for non edge-on galaxies. In a first study, we do not 
subtract the bulge in the image, and deproject it together 
with the disk, assuming that the bulge is as flattened as 
the disk. In a second study we subtract the bulge light 
from the image, before deprojecting the disk, and then we 
add it to the deprojected image, assuming the bulge is 
spherical. We believe that the two solutions will bracket 
the actual one, so we record the average value. However 
we checked that choosing one of the two solutions does not 
result in any significant change in our analysis of lopsid¬ 
edness: in the following we mainly study the disk surface 
density outside of the bulge radius, so that assumptions 
about the bulge component are not critical. 

A constant mass-to-light ratio and a constant scale- 
height are assumed. The value of m = 1 amplitude in 
surface density is computed for three values of the scale- 
height (225, 325, and 425 pc) and we keep the average 
value. 

2.3. The A\ parameter: Lopsidedness in surface density 
2.3.1. Definition 

The near-IR light is assumed to be the best tracer of 
old stars that dominate the stellar mass, less affected by 
dust extinction or star formation. We compute the Fourier 
transforms of the stellar surface density, taken from the 
NIR images. 

The surface density is decomposed as: 
p(r, cp) = po(r) + ^2 a m (r) cos (rrup - 4> m (r)) (1) 

m 

where the normalised strength of the Fourier component 
to is A m (r) = a m /po{r). Thus, A\ represents the normal¬ 
ized amplitude of the disk lopsidedness at a given radius, 
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Fig. 1. Radial variation of Ai, the normalized Fourier am¬ 
plitude of the stellar density, and phase 4>i ■, for the lopsided 
galaxy NGC 1637, measured from the near-infrared image 
of the OSUBGS sample. The disk scale-length is 2.45 kpc 
and the value of < A\ >= 0.19 was measured between 3.7 
kpc and 6.1 kpc. The near-IR image of the galaxy from the 
OSUBGS catalog is shown in Fig. M This galaxy shows 
a high degree of lopsidedness, but has no interacting com¬ 
panion and no evidence of a recent merger. 


and A 2 represents the same for an to = 2 distortion like 
a bar or spiral arms. The quantity 0 m denotes the posi¬ 
tion angle or the phase of the Fourier component to. We 
further define < A\ > as the average of the normalized 
amplitudes A\{r) between 1.5 and 2.5 exponential radii. 

Figure 1 shows a plot of the amplitude A±(r) and the 
phase 0i (r) versus radius r for NGC 1637, a typical lop¬ 
sided galaxy (Block et al. 1994); the deprojected density 
map computed from the near-IR image from the OSUBGS 
data is shown later in Fig. This illustrates the main 
features of disk lopsidedness, namely that A\ increases 
with radius, and the phase is fairly constant with radius 
in the outer parts. This is the case for most galaxies in 
our sample, even though there are a few exceptions, and 
it confirms the previous result by Rix & Zaritsky (1995) 
for a larger sample studied here. 

The distribution of A\ values vs. radius for the 163 
OSU sample of galaxies (with inclination < 70°) is given 



<A1> 


Fig. 2. Histogram of to = 1 Fourier amplitude in the stel¬ 
lar density deduced from the near-infrared images of the 
OSUBGS sample of galaxies. The < A\ > parameter is the 
average of the normalized Fourier amplitude for to = 1, 
over the radial range between 1.5 to 2.5 disk scalelengths. 


in the Appendix. The average values < A\ > and < A 2 > 
for the sample are listed in Table rm of Appendix El 

2.3.2. Statistical properties 

The histogram of the < A\ > parameter is displayed in 
Fig. 2, and the cumulative function 1 of < A\ > in Fig. 3. 
The histogram confirms and extends the previous results 
by Rix & Zaritsky (1995) for a larger sample. A good 
threshold for deciding that a galaxy is ’’lopsided” is < 
Ai >= 0.05: according to our numerical simulations in 
Sect. 4, a value of < A\ > larger than 0.05 is unlikely 
to result from internal mechanisms. Internal mechanisms 
may in theory trigger to = 1 asymmetries (e.g., Masset & 
Tagger 1997), but our simulations shown in Sect. 4, as well 
as other existing works, do not show strong lopsidedness 
arising spontaneously without any external perturbation. 
With this criterion, we find that 63% of the galaxies are 
lopsided. The mean value of < Ai > over our sample is 
0.11, and 34% of galaxies have < A\ > larger than this 
average value. 

In the above analysis, we have restricted the sample to 
galaxies inclined by less than <70°. This value is chosen 
because the statistical distribution of the < A\ > param¬ 
eter is nearly identical whether the cut-off is chosen to be 
30°, 50°, or 70° (see Fig. 3), and the higher value in this 
range was chosen as it gives a larger sample which is more 
useful for a statistical study. For higher values of the cut¬ 
off, statistical properties appeared to be changed, for the 
lopsidedness can no longer be correctly estimated. 


1 The cumulative function of < Ai > displays the fraction 
of galaxies with < Ai > larger than the x-axis value 
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Fig. 3. Cumulative functions of < A\ > (fraction of galax¬ 
ies with < Ai > larger than the x-axis value) for galax¬ 
ies seen below a certain cut-off in the inclination angle i\ 
for i < 30°, 50°, and 70°. The curves are nearly identi¬ 
cal, hence the highest value of the cut-off 70° is chosen, 
as it gives the largest number of galaxies in the sample. 
Galaxies more inclined than 70° no longer follow the same 
statistical distribution and have been rejected. 



Fig. 4. Histogram of m = 1 Fourier amplitude in the grav¬ 
itational potentials (Qi) for the OSUBGS sample of galax¬ 
ies. 

We prefer to use A i, the Fourier amplitude of lopsided¬ 
ness in surface density, as the main indicator of disk lop¬ 
sidedness since it is a directly measured quantity, and is 
frequently used in the literature. 


2.4. The Qi parameter 


We next derive the gravitational potential of each galaxy 
from the near-infrared images. This potential gives a 
global estimation of the perturbation of the galaxy disks, 
but is weighted quite differently from the density distri¬ 
bution. This gives a complementary view of the pertur¬ 
bations, weighted more globally, while the A\ parameter 
gives values which are weighted more locally. 

The image of the galaxy is completed in the 3rd di¬ 
mension, assuming an isothermal plane with a scale-height 
constant with radius and equal to 1 / 12th of the radial 
scale-length. The gravitational potential is deduced by the 
Fourier transform method on the grid assuming a constant 
mass-to-light ratio. 

The potential is decomposed as: 

$(r, 6) = $ 0 (r) + ^ ®™{r) cos (mO - <j) m ) (2) 


We define the strength of the m-Fourier component, 
Q m (r) = m$ m /r\Fo(r)\, and its global strength over the 
disk as max,. (Q m (r)) (e.g., Combes & Sanders 1981). 
The maximal torque at a given radius is defined by 


Qt{R) = 


jpmax 


(R) 


F 0 (R) 


1 / d$(R,6) \ 

R v dO ) max 

d® 0 (R) 

dR 


(3) 


where F™ ax (R) represents the maximum amplitude of the 
tangential force at radius R 1 and F 0 (R) is the mean ax- 
isymmetric radial force inferred from the m =0 component 
of the gravitational potential. Our interest is focussed here 
on the m = 1 component, Q\, defined as the maximum of 
Qi(R) between 1.5 and 2.5 disk scale-lengths. 

We display the histogram of the values of Q\ in Fig. 4. 
This is similar to the histogram for < A\ > as expected. 


2.5. Lopsidedness of gaseous disks 

The asymmetries in the HI component have been esti¬ 
mated in a large sample on global velocity profiles by 
Richter & Sancisi (1994). Since galaxies with interfero¬ 
metric HI maps are rare, they investigate the usability of 
global HI profiles as tracers of the lopsidedness of the HI 
disks in galaxies and notice that asymmetries in the global 
profiles usually correlate with those found on the HI maps. 
There are, however, cases where the HI map asymmetry 
does not translate into the HI global profile asymmetry, 
because of projection effects of the velocity along the line- 
of-sight, or because the asymmetry changes sign over the 
disk and conspires to reduce the global projected the per¬ 
turbation. Therefore the velocity profile tracer is likely to 
underestimate the asymmetry level. 

We gathered HI velocity profiles from the literature for 
our present sample, in order to estimate the neutral gas 
lopsidedness. The HI asymmetry was measured in terms 
of the second criterion from Richter & Sancisi (1994). We 
first compute the flux ratio / between the two halves of 
the profiles (the two horns), as illustrated in Fig. 0 to 
compare the mass of the approaching side of the galaxy, 
and the mass of the receding side. To measure the flux 
of each horn, we took the total width at 20% of peak flux 
density, and determined both the extremities of the profile 
and its center, i.e. the division between the two halves. 

We define the asymmetry parameter as = 10(1 — /), 
where / is the flux ratio of the two horns of the profile. 
HI data available for 76 galaxies of our sample. The corre¬ 
sponding histogram of E\ is shown in Fig. EO which shows 
that asymmetries in the HI component are very frequent. 
The asymmetry is larger than 10% (E\ larger than 1) 
for nearly two thirds of the galaxies in our sample, while 
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Velocity (km/s) 

Fig. 5. HI velocity profile of NGC 150, showing the two 
horns, the center of the profile and the extremity of the 
two horns. Fi and F 2 are the integrated flux of each horn, 
and their ratio gives an estimate of the HI asymmetry. 



El 


Fig. 6. Histogram of the E\ parameter measured on HI 
velocity profiles. 

in the stellar component only one third of galaxies have 
< Ai > larger than 0.1. Asymmetries of more than 20% 
(< A\ >> 0.2 or Ei > 2) are also much more common 
in the HI component than in NIR disks. Even if E\ only 
provides a lower limit to the actual HI asymmetry, we find 
that in our sample lopsidedness is even more frequent and 
stronger in gas than in stars. 

3. Dependence of the lopsidedness on various 
physical parameters 

In this section, we study the statistical relations between 
lopsidedness and other parameters (the Hubble type, the 
presence of companions, and the presence of m = 2 asym¬ 
metries like arms or bars). The results presented here are 
for the sample of galaxies below the inclination cut-off of 
70°, but we made sure that the results are not artifacts 
of deprojection by checking that they remain unchanged 
for a cut-off of 50° and even for a cut-off of 30° when the 
number of galaxies below 30° was large enough to make 
this study possible. 



Fig. 7. Cumulative function of < A\ >, for three groups 
of Hubble types of spiral galaxies: the early-types (0 < 
T < 2.5), the intermediate-types (2.5 < T < 5), and the 
late-types (5 < T < 7.5). The late-type spirals are more 
lopsided than the early-type ones. 

3.1. Dependence of < A\ > on galaxy type 

To study the dependence of the lopsidedness < A± > on 
galaxy types, we plot the cumulative functions of < A\ > 
for three groups of spiral galaxies in Fig. 7 : the early-types 
(0 < T < 2.5), the intermediate-types (2.5 < T < 5), and 
the late-types (5 < T < 7.5). The correspondence between 
the Hubble type parameter T and the Sa-Sd classification 
of early-type/late-type spiral galaxies, as explained for in¬ 
stance by Binney & Merrifield (1998), is: T ~ 2 for Sa, 
T ~ 4 for Sb, T ~ 7 for Sc. The number of galaxies with 
T > 7.5 in our sample is not large enough to give any 
representative result on the extremely late-type systems. 

The late-type galaxies have a higher fraction of lop¬ 
sided galaxies, as can be seen for any value of lopsided¬ 
ness < Ai >. This result is similar to what Matthews et 
al. (1998) found in HI data. 

If we had employed Q 1 instead of Ai, we could think 
that this correlation between the asymmetry and the 
Hubble type is an artefact of bulge contamination, since 
the axisymmetrical contribution of the bulge to the poten¬ 
tial reduces the value of Qi in early-type spiral galaxies 
more than in late-type galaxies, even beyond the bulge ra¬ 
dius. However, since we have used the Ai parameter and 
taken its average value between 1.5 and 2.5 disk scale- 
length, our results are not affected by bulge contamina¬ 
tion, except maybe in a few cases were the bulge has a 
large radial extent, but this is rare. Then, the correlation 
between < Ai > and the Hubble type really means that 
late-type spirals are more lopsided than early-type ones. 

3.2. Comparison of < A\ > and < A 2 > 

In Fig. 8, we show the cumulative distribution of < A\ > 
for galaxies that have strong and weak m = 2 asymme¬ 
tries, and find that there is a clear correlation between 
Ai and A 2 . A similar result is obtained when we mea- 
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Fig. 8. Cumulative function of < A\ > for galaxies with 
values of < H 2 > higher than the median value of < A 2 > 
over our sample (dashed line), and < A 2 > smaller than 
this threshold (solid line). Galaxies with large m = 2 
asymmetries are clearly more lopsided than galaxies with 
small in = 2 asymmetries. Here < A 2 > is the average 
value of A 2 between 1.5 and 2.5 disk scalelengths. The 
correlation between A\ and A 2 remains if we take the av¬ 
erage value of A 2 between 0.5 and 1.5 scalelengths (which 
would rather correspond to a central bar), or between 2.5 
and 3.5 scalelengths (which would rather correspond to 
outer spiral arms). 

sure < A 2 > as the average value between 0.5 and 1.5 
disk scalelengths or between 2.5 and 3.5 disk scalelengths, 
suggesting that the presence/strength of a lopsidedness is 
correlated to the presence/strength of m = 2 asymme¬ 
tries, whatever the exact nature of this m = 2 asymmetry 
(central bar and or outer spiral arms). 

3.3. Dependence of A\ on the presence of companions 


0.4 

A 

<0.3 

v 

0.2 


0.1 -1-1-1-1-1- 

-7 -5 -3 -1 

■p 

Fig. 9. Correlation between the < A\ > parameter and 
the tidal parameter Tp defined in text, for the 35 most 
lopsided galaxies with i < 70°. No correlation is found. 
Strong lopsidedness in very isolated galaxies (Tp < —6) 
is frequent, and is not the privilege of galaxies with com¬ 
panions. 

The absence of correlation between strongly lopsided disks 
and the presence of companions agrees with what Wilcots 
& Prescott (2004) found. Note that at this stage, it does 
not rule out interactions/mergers as being responsible for 
disk lopsidedness: the explanation could be that lopsided¬ 
ness is very long-lived. 

4. Origin of disk lopsidedness 

To understand the origin of lopsided disks, what must be 
explained is not only the frequency and strength of lop¬ 
sidedness, but also the dependence on the Hubble type 
and m = 2 asymmetries and the non-correlation of strong 
lopsidedness with the presence of interacting companions. 


To quantify the effect of tidal forces, we define the tidal 
parameter as 



where the sum is computed over the companions, Mq is 
the mass of the galaxy considered, and i?o its scalelength. 
Mi is the mass of each companion, and Di its projected 
distance on the sky. The ratio of masses is estimated from 
the ratio of blue luminosities. The search of neighbors was 
done with the NED database, with a search radius of 2.5 
degrees on the sky. We take only the companions with 
measured radial velocities within 500 km s . 

The < Ai > parameter values for the most lopsided 
galaxies are displayed versus the tidal parameter in Fig. 9. 
No correlation is found. Strong lopsidedness with < A\ > 
as large as 0.1 or larger is frequent in isolated galaxies 
(with Tp < —4) without any close companion or any sign 
of a recent merger. NGC 1637 is a well-known example, 
but there are many other similar objects in our sample. 


4.1. Scenarios for the origin of lopsidedness 

Minor m = 1 perturbations can spontaneously arise in 
a disk. For instance, a two-arm spiral structure is never 
perfectly symmetrical, which results in a non-zero value of 
< Ai >. The numerical simulations detailed below, when 
run for isolated galaxies, show that the value of < A\ > 
arising spontaneously in an isolated galaxy rarely exceeds 
0.05 (in simulations of isolated galaxies, < Ai > is smaller 
than this threshold 88% of the time), and is never as large 
as 0.1 between 1.5 and 2.5 disk scale-lengths (the inner 
lopsidedness, at small radii, can be stronger, but is not 
studied in this paper). Thus, the lopsidedness observed 
in most spiral galaxies is the result of an external per¬ 
turbation. Three scenarios to explain a galaxys becoming 
significantly lopsided can then be proposed: 

— distant interaction with a galaxy of comparable mass. 
The companion must have a high velocity and/or be 
distant enough for this encounter not to result in a 
merger; otherwise the merger remnant would be an 
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Type 

bulge-to-disk 
mass ratio 

halo-to-disk 
mass ratio 

Sa 

0.5 

0.4 

Sb 

0.3 

0.5 

Sc 

0.2 

0.6 

Sd 

0.1 

0.8 


Table 1. Bulge, disk, and halo parameters used in nu¬ 
merical simulations. The halo mass that is accounted for 
is the mass of dark matter inside the disk radius. 


elliptical galaxy, while we only have disk galaxies in 
our sample. 

— a minor merger with a small companion, with a mass 
ratio such as 10:1 or even more: such a companion is 
not massive enough to transform the system into an 
elliptical galaxy, but may induce important perturba¬ 
tions like lopsidedness. 

— accretion of external gas. If the accretion is asymmet¬ 
rical, this may result in a lopsided disk, even in the 
stellar component, when the gas is converted into stars. 

We now study these three scenarios in numerical simu¬ 
lations, and compare the results of each scenario with the 
properties of the observed galaxy sample. 

4.2. Numerical techniques 

In our N-body code, galaxies are described by particles 
of stars, gas, and dark-matter. We use N = 10 6 par¬ 
ticles. The gravitational interactions are computed with 
the three-dimensional FFT code of Bournaud & Combes 
(2003) on 256 3 grids, with a spatial gravitational softening 
of 300 pc. The dissipative nature of the ISM is modeled by 
the sticky-particles scheme of Bournaud & Combes (2002), 
with elasticity parameters (3 r = fit = 0.8. Star formation 
is described by a generalized Schmidt law (Schmidt 1959), 
and time-dependent stellar mass-loss is described by a 
model inspired by Jungwiert et al. (2001). More details 
about these codes can be found in Bournaud & Combes 
(2002 and 2003). 

Each spiral galaxy is modeled by a stellar and gaseous 
Toomre disk of scale-length 5 kpc and radius 15 kpc, of 
mass 10 11 Mg, a central bulge, and a spherical dark halo. 
All the components are live (treated with self-gravitating 
particles) in these simulations. The bulge-to-disk and 
dark-to-visible mass ratios that we have used, as a func¬ 
tion of the Hubble type (Sa-Sd), are given in Table|T] The 
run parameters for each simulation are given later. 

In these simulations we measure < A\ >, < A -2 >, Qi, 
and Q 2 , at each time step, exactly as in observed systems, 
so that a robust comparison can be made. 

4.3. Distant galactic encounters 

We first study distant encounters (without mergers) be¬ 
tween galaxies of comparable masses. The parameters of 
the simulations that we ran are shown in Table □ We 


focused on situations favorable to the formation of large 
lopsidedness in the disk. This is mainly the case on retro¬ 
grade orbits, which is naturally explained by the fact that 
the pattern speed of m = 1 asymmetries, fl — k, is neg¬ 
ative in a galactic disk, in contrast to to = 2 spiral arms 
or bars that are more easily triggered on direct orbits, for 
fl — k/2 is positive (e.g., Gerin et al. 1990). 



Time (Gyr) 


Fig. 10. Simulation of distant interaction between two 
galaxies of comparable mass (Run 9). Top: temporal evo¬ 
lution of < Ai > (solid line) and the tidal parameter Tp 
(dashed, defined in Sect. 3). A strong lopsidedness is trig¬ 
gered during the interaction. It is observed for about 4 
Gyrs, but is much lower once the system is isolated. - 
Bottom: Q\ parameter in the disk plane for the total 
potential (solid line), the internal target galaxy poten¬ 
tial only (dashed curve), and the external potential of the 
companion only (dotted curve). During the interaction, a 
strong lopsidedness results from the potential of the com¬ 
panion, but is not totally supported by the potential of 
the target galaxy itself. The maximal lopsidedness in Q 1 
can be decomposed as an external contribution from the 
companion (short-lived) and an internal lopsidedness sup¬ 
ported by the target galaxy potential (long-lived). This 
explains why there is a strong peak in Q\ and A\ during 
the interaction, and a long-lived but weaker lopsidedness 
after the interaction. 

The temporal evolution of the < A\ > parameter for 
a representative case (Run 9) is shown in Fig. HOI During 
the interaction itself, a strong lopsidedness with values of 
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N° 

M 

T 

Orient. 

V (kin s 1 ) 

r (kpc) 

e 

1 

1:1 

Sb 

D 

200 

220 

0 

2 

1:1 

Sb 

D 

300 

130 

0 

3 

1:1 

Sb 

D 

150 

350 

0 

4a 

1:1 

Sa 

R 

250 

130 

0 

4b 

1:1 

Sb 

R 

250 

130 

0 

4c 

1:1 

Sc 

R 

250 

130 

0 

5a 

1:1 

Sa 

R 

160 

300 

0 

5b 

1:1 

Sb 

R 

160 

300 

0 

5c 

1:1 

Sc 

R 

160 

300 

0 

6a 

1:1 

Sa 

R 

200 

220 

0 

6b 

1:1 

Sb 

R 

200 

220 

0 

6c 

1:1 

Sc 

R 

200 

220 

0 

7 

1:1 

Sb 

R 

250 

130 

60 

8 

1:2 

Sb 

R 

250 

130 

0 

9 

2:1 

Sb 

R 

250 

130 

0 

10 

1:3 

Sb 

R 

250 

130 

0 

11 

3:1 

Sb 

R 

250 

130 

0 

12 

1:2 

Sb 

R 

160 

300 

0 

13 

2:1 

Sb 

R 

160 

300 

0 

14 

1:2 

Sb 

R 

450 

450 

0 

15 

2:1 

Sb 

R 

450 

450 

0 


Table 2. Parameters for the simulations of distant inter¬ 
actions of galaxies. M is the mass ratio (1:2 means that 
the companion is smaller than the target galaxy, 2:1 that 
it is more massive). T is the Hubble Type (see Tabic EJ. 
The orientation of the orbit can be direct (D) or retro¬ 
grade (R). We also give the relative velocity at an infinite 
distance V, the impact parameter r, and the inclination 
of the orbital plane with respect to the galactic plane 9. 


< Ai > higher than 0.15 is triggered. It is short-lived, 
lasting for about 500 Myrs. After that, a moderate but 
long-lived m = 1 mode is observed during at least 4 Gyrs. 
All the simulations of distant galaxy interactions show the 



Fig. 11. Largest value of < A\ > as a function of the 
tidal parameter Tp for all the simulations of tidal galaxy 
interactions. All the results of the simulations lie below 
the solid line in this (Ai;Tp) plane. This scenario does 
not explain the presence of strong lopsidedness (< A\ > 
larger than 0.1) in isolated galaxies (tidal parameter < 
—5), while there are many observed cases. 


same behavior. This result is explained by analysing of the 
potential in terms of the Q i parameter, shown in Fig. m 
The strong peak in < A\ > and Q \ observed during the 
interaction is supported for one part by the internal po¬ 
tential of the target galaxy, and for the other part by the 
external potential of the companion that pulls-up the disk. 
When the companion has left, only the part supported by 
the target galaxy potential remains; this mode is long- 
lived, but weaker than the peak observed in Q\ and in 
Ai, during the interaction. 

As a consequence of this temporal evolution, a strong 
lopsidedness is only observed when the companion is still 
close to the galaxy. This is the case in all the simulations, 
as shown in Fig. El with the largest value of < A\ > over 
our whole simulation sample as a function of the tidal pa¬ 
rameter. Even over the whole parameter space that we 
have covered, we never find a lopsidedness that is long- 
lived enough and strong enough to result in < A\ > val¬ 
ues higher than 0.1 in galaxies that appear very isolated. 
Strong lopsidedness, with < A\ > of 0.15 or more, are 
found only when the tidal parameter Tp is larger than —4 
in this scenario: depending on the mass ratio, this typi¬ 
cally implies that the companion is within a few hundreds 
of kpc from the galaxy, and has not left the scene yet. 

4.4. Minor galaxy mergers 

The second scenario that we study is the merger with a mi¬ 
nor companion. Walker et al. (1996) have already found 
that this process can form significant disk lopsidedness, 
visible for at least 1 Gyr in their simulations . The corre¬ 
sponding run parameters are given in Table |oJ Again, the 
most favorable cases are on retrograde orbits. A mass ratio 
of 10:1 can trigger significant lopsidedness. Mergers with 
smaller mass ratios, such as 7:1, are even more efficient 
for triggering m = 1 modes, but they will also lead to ma¬ 
jor disturbances of the system (see Bournaud et al. 2004, 
2005), making the galaxy become an extremely early-type 
spiral or even an SO-like object, so they are not good can¬ 
didates for the formation of lopsidedness in spiral galaxies 
in general, and we do not simulate mass ratios larger than 
7:1. At the opposite, large mass ratios like 20:1 trigger 
weaker lopsidedness with typical values of < A\ > smaller 
than 0.1. We have not studied mass ratios larger than 20:1 
that cannot be the clue to the frequent lopsidedness with 
< Ai » 0.1 that is observed. 

The temporal evolution of < Ai > in a minor merger 
(see Fig. El for runs 4b and 6c) is fairly similar to what 
occurs in distant galactic encounters without mergers. A 
very strong lopsidedness is triggered during the merger 
with values of < A\ > that can be larger than 0.2, but 
this lasts typically no more than 500 Myrs. Later on, a 
long-lived m = 1 mode is still present, but is weaker with 
values of < Ai > not much larger than 0.1. The strong 
peak in < Ai > observed during the merger is mainly 
related to the very asymmetric distribution of the stars 
from the companion; we find a lower value of < Ai > when 
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Fig. 12. Temporal evolution of < A\ > in two minor 
mergers (top: Run 4b, bottom: Run 6c). The solid line 
is when stars from the companion (or formed in gas from 
the companion) are included. The dashed line corresponds 
only to stars from the target galaxy (or formed in gas 
clouds present in the target before the merger). When the 
two curves have significantly different values, this means 
that the system is highly disturbed, for instance at t\ (see 
the corresponding snapshot for Run 4b on Fig. IIdll . Tidal 
tails, streams, and other merger indicators, are still seen 
at f 2 - The system can be classified as an ’’isolated galaxy” 
after the two curves have reached similar values, for in¬ 
stance at to (see Fig. 113. We then consider that the max¬ 
imal lopsidedness for an ’’isolated” galaxy is found when 
the two curves reach similar values, at to' corresponding 
values of < A\ > are given for each simulation in Table |3 


we only consider the target galaxy material (see Fig. m- 
After 2-3 rotations, the companion stars are mixed in the 
system and follow the same azimuthal distribution as the 
other stars - even if their radial distribution is different. 
The value of < Ai > is then smaller, but present in all 
the stars, regardless of whether stars from the companion 
are included in the analysis or not. 

Can minor mergers create a strong lopsidedness when 
the galaxy appears ’’isolated”, once the most obvious 
merger-related features have disappeared? A strong lop¬ 
sidedness does arise only during the merger, but the 


N° 

M 

T 

Orient. 

V 

km s -1 

r 

kpc 

e 

< Al > 
isolated 

1 

10:1 

Sb 

D 

100 

50 

0 

0.08 

2a 

10:1 

Sa 

R 

100 

50 

0 

0.11 

2b 

10:1 

Sb 

R 

100 

50 

0 

0.12 

2c 

10:1 

Sc 

R 

100 

50 

0 

0.11 

3 

10:1 

Sb 

D 

150 

30 

0 

0.07 

4a 

10:1 

Sa 

R 

150 

30 

0 

0.11 

4b 

10:1 

Sb 

R 

150 

30 

0 

0.10 

4c 

10:1 

Sc 

R 

150 

30 

0 

0.11 

5 

10:1 

Sb 

D 

60 

40 

0 

0.09 

6a 

10:1 

Sa 

R 

60 

40 

0 

0.11 

6b 

10:1 

Sb 

R 

60 

40 

0 

0.13 

6c 

10:1 

Sc 

R 

60 

40 

0 

0.12 

7 

20:1 

Sb 

D 

100 

50 

0 

0.07 

8a 

20:1 

Sa 

R 

100 

50 

0 

0.10 

8b 

20:1 

Sb 

R 

100 

50 

0 

0.10 

8c 

20:1 

Sc 

R 

100 

50 

0 

0.12 

9 

20:1 

Sb 

D 

150 

30 

0 

0.08 

10 

20:1 

Sb 

R 

150 

30 

0 

0.11 

11 

10:1 

Sb 

R 

60 

40 

60 

0.09 

12 

10:1 

Sb 

R 

100 

50 

60 

0.10 

13 

20:1 

Sb 

R 

200 

30 

0 

0.09 

14 

10:1 

Sb 

R 

200 

50 

0 

0.11 

15 

7:1 

Sb 

R 

60 

40 

0 

0.13 

16 

7:1 

Sb 

R 

100 

50 

0 

0.13 

17 

7:1 

Sb 

R 

150 

30 

0 

0.14 


Table 3. Parameters for the simulations of minor mergers. 
M is the mass ratio. T is the Hubble Type, determining 
several bulge/halo parameters (see Table ^ f° r details). 
The orientation of the orbit can be direct (D) or retro¬ 
grade (R). We also give the relative velocity at an infinite 
distance V , the impact parameter r, and the inclination 
of the orbital plane on the galactic plane 9. The last col¬ 
umn indicates the value of < A\ > when the stars of the 
companion are azimuthally mixed in the galaxy (i.e. when 
the value of < A\ > is unchanged when we include these 
stars or not, at the instant called to on Fig. [Hi! . Before 
that, the lopsidedness can be stronger, but the system ap¬ 
pears like an ongoing merger, not like an isolated galaxy, 
as explained in text and shown in Figs. d and El The 
value indicated here corresponds to the instant at which 
the galaxy begins to appear ” isolated”. 


galaxy is then largely disturbed and would be classified as 
an ” on-going merger”, rather than as an ” isolated galaxy”. 
The question is then to know whether a strong lopsided¬ 
ness can still be present when the merger is relaxed and the 
galaxy looks ’’isolated”. As explained, the highest values 
of < Ai > are observed when the azimuthal distribution 
of the stars of the companion is very different from the 
distribution of the other stars, at the instant called ti on 
Fig. El This big difference in the two stellar components 
is enough to attest that a merger is occurring, and the 
visual aspect of the system clearly shows that it is an on¬ 
going merger, not an isolated object (see Fig. HI- Even at 
time t- 2 , when including or not the stars from the compan¬ 
ion does not largely change the value of < A\ >, there are 
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Fig. 13. Snapshots of the minor merger scenario, Run 4b (see the temporal evolution of < A\ > and the definition of 
ti, < 2 , t%, in Fig. El- At t\, there is a strong lopsidedness, but the ongoing merger is clearly visible. Even at £ 2 , when 
the stars from the companion are more ’’mixed” in the system, and no longer dominate the value of < A\ >, there 
are still obvious signs of the recent merger (like arcs and tidal debris in the outer parts, that are even more visible in 
the gas component). Only at < 3 , the system is an ’’isolated” galaxy without sign of a past merger, but the amplitude 
of the lopsidedness is much lower, as shown by Fig. IT21 


still indicators that a recent merger has occurred, such as 
streams or tidal debris (see Fig. m- It’s only at time £ 3 , 
a few dynamical times later, that the system looks like an 
” isolated” galaxy - but the value of < A\ > is then twice 
lower than at t\. Thus, the galaxy can look ’’isolated” only 
after we find similar values of < A\ > with and without 
including the stars from the companion. Before this time, 
there is evidence that it is an ” ongoing merger”. Then, in 
Table |21 we included in the value of < A\ > at the instant 
when the stars from the companion are azimuthally mixed 
in the system (called f 0 in Fig. 1121 . The result is that, 
while the value of < Ai > can be larger than 0.20 during 
the merger itself, the lopsidedness is much more modest 
once the merger remnant looks relaxed and can be clas¬ 
sified as an ’’isolated” galaxy: at this moment, the value 
of < Ai > is typically ~ 0.10 and is always smaller than 
0.15, even for the not-so-minor 7:1 mergers - mergers with 
companions more massive than 7:1 are not considered, for 
the merger remnant is not a typical spiral-like system any 
more (see Bournaud et al. 2005). 

4.5. Difficulties with the galaxy interaction/merger 
scenarios 

Galaxy interactions, either distant encounters or minor 
mergers, can trigger strong lopsidedness. The galaxy can 
still be lopsided a few Gyrs after the event. However, sev¬ 
eral observational results cannot be explained by these 
scenarios. 

Isolated galaxies that are strongly lopsided. As mentioned 
in Sect. 3, no correlation has been found between the 
strongest lopsidedness and the presence of interacting 
companions. Strong lopsidedness with < A\ >> 0.1 are 
frequently found in galaxies that are very isolated and do 
not show any sign of a recent merger. 


This does not directly rule out galaxy interactions as 
responsible for lopsided disks, because a strong lopsided¬ 
ness may a priori still be observed long after the interac¬ 
tion/merger, when no sign of it is observable any more. 
However, our simulations show that in the case of dis¬ 
tant interactions with massive companions, the strongest 
lopsidedness can only be observed during the interaction, 
when the companion is still rather close, so that the strong 
lopsidedness with < A\ >> 0.1 in very isolated galaxies 
(tidal parameter of —8 to — 6 ) that are frequently ob¬ 
served cannot be explained; on the contrary, with this 
scenario, a correlation of < A\ > with the tidal param¬ 
eter is expected, at least for the most lopsided galaxies. 
Minor mergers do not provide a better explanation, since 
the lopsidedness is very strong only during the merger it¬ 
self, when the stellar population of the companion is not 
yet mixed in the disk, with other visible indicators that a 
merger has occurred recently (tidal tails, streams from the 
companions, etc.). Strong lopsidedness in galaxies that do 
not show any evidence of a recent merger and no massive 
interacting companion (like NGC 1637, M 101, and many 
others) cannot be accounted for. 

Correlation with m = 2 asymmetries. We have shown in 
Sect. 3 that lopsidedness is correlated with the strength 
of m = 2 asymmetries (spiral arms and/or bars). This 
result could have no relation with the physical origin of 
lopsidedness, for it could simply be that both lopsided¬ 
ness and m = 2 asymmetries are spontaneously stronger 
in the same type of galaxies. However there are at least 
two pieces of evidence that it is not the case: (z) lopsided¬ 
ness is stronger in late-type galaxies (Sect. 3 and Fig. 7), 
while m = 2 spiral arms/bars are stronger in early-type 
spiral galaxies (e.g. Elmegreen & Elmegreen 1985, and also 
our own sample, see Fig. m, and (zz) when we take only 
early-type or late-type spirals, the dependence between 
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Fig. 14. Cumulative function of < A 2 > for late-type, 
intermediate-type, and early-type spiral galaxies (a similar 
figure for < A\ > has been shown in Fig. 0 . Here < 
A 2 > is the average value of A 2 between 1.5 and 2.5 disk 
scalelengths, but the anti-correlation of < A 2 > with T 
remains if we consider the average values between 0.5 and 
1.5, or 2.5 and 3.5 disk scalelengths. Our sample obeys the 
general rule that m = 2 bars/arms are stronger in early- 
type spiral galaxies. This is the opposite of what we found 
for lopsidedness, which is stronger in late-type galaxies. 

A\ and A 2 remains as strong as for the whole sample (see 
Fig. II 5ll . The correlation between Ai and A 2 is thus re¬ 
lated to the formation of lopsidedness and m = 2 asym¬ 
metries. Yet, galaxy interactions/mergers do not explain 
it because: (?) the strongest lopsidedness are triggered on 
retrograde orbits, while (??) m = 2 asymmetries are trig¬ 
gered on prograde orbits; on retrograde orbits they are 
not triggered but can instead be weakened (e.g., Gerin et 
al. 1990, Noguchi 1996). One would then expect weaker 
m = 2 asymmetries in the most lopsided galaxies, which 
contradicts of the observational reality. The dependence 
observed between A\ and A 2 cannot be explained if one 
assumes that galaxy interactions are responsible most of 
the observed lopsidedness. 

Correlation with Hubble type. We have shown that there 
is a dependence of A\ on the galaxy type: late-type spiral 
galaxies are much more lopsided than early-type ones. The 
same result has been found in gaseous disks by Matthews 
et al. (1998). This could mean that it is easier to trigger a 
strong lopsidedness in a late-type galaxy. However, simula¬ 
tions of interactions and mergers do not show significantly 
larger lopsidedness in late-types (there might be a small 
dependence on the Hubble type in simulations, but not 
as large as the observed one). Then, lopsidedness does not 
seem to be intrinsically stronger in late-type galaxies. The 
dependence on the observed Hubble Type must be a result 
of the formation mechanism, but it cannot be explained 
by galaxy interaction scenarios. 

Lopsidedness is not only stronger in late-type galax¬ 
ies, where there are more values of < Ai > ~ 0.15 — 0.20 



Fig. 15. Same as Fig. [S] for early-type spiral galaxies only 
(left, 0 < T < 4) and late-type galaxies only (right, 4 < 
T < 8 ). As discussed in the text, this shows that the 
correlation between m = 1 and m = 2 asymmetries found 
for the whole sample in Sect. 3 does not just mean that 
they are stronger in the same types of galaxies, but really 
means that they are preferentially formed simultaneously. 
Indeed, this correlation is still found when we limit the 
sample to early- or late-type galaxies. 


than in early-type (see Fig. 7), they are also more fre¬ 
quent in late-type galaxies: if we choose the threshold for 
a galaxy to be regarded as lopsided at < A\ >= 0.05, 
as suggested by our simulations, then 96% of late-type 
galaxies are lopsided, while only 62% of early-type spi¬ 
ral disks are lopsided (and this difference still holds if we 
chose another threshold). While galaxy interactions could 
have a priori explained that lopsidedness is stronger in 
late types (even when the simulations finally show they 
do not), they cannot explain at all why it is more fre¬ 
quent. On the contrary, a galaxy interaction makes the 
system become earlier-type. For the minor mergers, that 
are much more likely than distant equal-mass interactions, 
we have shown that the merger must not be too minor to 
trigger strong lopsidedness: mass ratios such as 10:1 or 
7:1 are efficient, but 20:1 mergers are already much less 
efficient. According to our recent study of galaxy mergers 
(Bournaud et al. 2005), a 10:1 merger is enough to make 























12 


Bournaud et al.: Lopsided spiral galaxies: evidence for gas accretion 


the system evolve significantly towards early-type (more 
massive bulge, thicker disk, less gas), especially on retro¬ 
grade orbits; and a 7:1 merger can even convert a typical 
spiral galaxy into a SO-like object. This evolution towards 
early-type is significant enough to be observable. Indeed, 
the variations in the bulge masses, velocity dispersions, 
and disk thickness, that result from a 10:1 merger corre¬ 
spond to a decrease in the Hubble type indicator T by a 
few units. So, minor mergers would at the same time trig¬ 
ger an m = 1 asymmetry and make the system become 
earlier-type. Distant galaxy interactions also would trig¬ 
ger lopsidedness and make the system become significantly 
earlier-type at the same time. 

Thus, the galaxy interaction/merger scenarios would 
tend to create more lopsidedness in early-type galaxies, 
because even if the galaxy is not early-type before the in¬ 
teraction, it will become early-type during the interaction 
itself, especially for the strongest lopsidedness for which a 
strong interaction/merger is required. An anti-correlation 
between < A\ > and the Hubble type T would then be 
expected if all the disk lopsidedness were the results of 
galaxy interactions/mergers. But the observational result 
is that lopsidedness is stronger and more frequent in late- 
type galaxies, which the interaction/merger scenarios fail 
to explain. 

Conclusion Several properties of lopsided galaxies, derived 
from our observed sample in Sect. 3, cannot be explained 
if one assumes that disk lopsidedness is mostly the result 
of galaxy interactions or mergers. 

We do not deny that distant galaxy interactions and 
minor mergers are responsible for the observed lopsided¬ 
ness in some galaxies. This is shown by our simulations 
as well as previous works; and since interactions/mergers 
do occur in the Universe, they have triggered some of the 
observed lopsidedness. But we find that all the lopsided 
disks cannot be the result of these processes, for otherwise 
many observational properties could not be explained at 
all. We cannot give the exact fraction of lopsidedness that 
is not the result of galactic interactions, but it must be 
large enough to account for the overall statistical proper¬ 
ties of our sample, so it is certainly several tens of percent. 

4.6. Gas accretion 

Simulations of gas accretion Since a large fraction of lop¬ 
sided disks are not the result of galaxy interactions and 
mergers and not the results of internal processes (that 
rarely lead to < Ai > values larger than 0.05, as men¬ 
tioned in Sect. 4.1), another scenario for making a galaxy 
lopsided must be proposed. We can then imagine that lop¬ 
sidedness is the result of gas accretion. Standard cosmolog¬ 
ical models predict that galaxies accrete large amounts of 
gas along cosmological filaments (e.g., Keres et al. 2004). 
The mass accreted over a few billion years represents a 
significant fraction of the disk mass itself. This accretion 
not only fuels the gaseous disk, but also the stellar one, 



Fig. 16. Configuration of accretion simulations with one 
filament. The spatial offset of the filament with respect to 
the disk center is given in TableQ]for each simulation. The 
size of the filament is 25 kpc. 


N° 

T 

M 

offset 

mean A\ 

1 

Sb 

6 

15 kpc 

0.12 

2 

Sb 

6 

25 kpc 

0.16 

3 

Sb 

6 

40 kpc 

0.15 

4 

Sb 

3 

25 kpc 

0.11 

5a 

Sa 

9 

25 kpc 

0.18 

5c 

Sc 

9 

25 kpc 

0.21 


Table 4. Parameters for simulations of accretion: initial 
Hubble type of the galaxy, accretion rate (in Mg yr -1 ), 
and spatial offset of the filament with respect to the disk 
center. The last column indicates the average value of < 
A\ > between t = 2 and t = 4 Gyrs. 


through star formation. It has already been shown that 
gas accretion fueling both the gaseous and stellar disk can 
be responsible both for the maintenance of spiral arms 
during a Hubble time, while spiral arms in isolated mod¬ 
els generally do not last more than 2-3 Gyrs, and for the 
reformation of bars if they have been dissolved either by a 
merger or by an internal mechanism (Bournaud & Combes 
2002, Block et al. 2002). 

This accretion can be asymmetrical, since distributed 
along a few filaments, and the accretion rate can be dif¬ 
ferent for each filament. This can then trigger an m = 1 
asymmetry in the disk. In this paper we do not try to 
simulate ’’realistic” accretion conditions predicted in de¬ 
tail by any cosmological model, but simply to understand 
whether asymmetrical accretion of gas and subsequent 
star formation can create a significantly lopsided disk. To 
this aim, we first consider accretion along one single fila¬ 
ment, off-centered with respect to the galaxy (see Fig.Hol). 
and will discuss more realistic cases with several filaments 
next. The velocity distribution of the accreted gas is as¬ 
sumed to be the same as in the outer disk. We ran simu¬ 
lations for accretion rates of 3, 6, and 9 Mg-yr -1 . Typical 
values of < A\ > in this model are given in Table 01 The 
temporal evolution of < A\ > is shown in Fig. ^3 for two 
cases. 

The main result of these simulations is thus that asym¬ 
metrical accretion of gas can create strong m = 1 asym¬ 
metries, even in the stellar component after the formation 
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Fig. 17. Top: temporal evolution of A\ in a simulation 
of accretion (Run 2). A strong lopsidedness is triggered, 
even in the stellar component (through star formation), 
in a galaxy that looks completely isolated. Then the value 
of < Ai > decreases slowly while the disk re-orients 
itself, but remains larger than 0.1 during more than 3 
Gyrs. Bottom: same simulations with no accretion af¬ 
ter t = 2 Gyr. The lopsidedness is long-lived (about 3 
Gyrs), but slightly shorter-lived than in the galaxy inter¬ 
action/merger scenario, probably because perturbations 
to the halo are weaker in this case. 

of stars. The degree of lopsidedness decreases after a few 
Gyrs, even if accretion is not stopped, because the disk re¬ 
orients itself in regard to the filament. When accretion is 
very asymmetrical, as in our model with one filament, the 
values of < A\ > can be as strong as in the extreme ob¬ 
served cases. This single-filament model is very arbitrary, 
and for a real galaxy it is more likely that there are sev¬ 
eral filaments with different accretion rates (e.g., Semelin 
& Combes 2005). Then, the lopsidedness will generally be 
more moderate; however, it will remain true in such a re¬ 
alistic situation that asymmetrical gas accretion results 
in a lopsided disk. Indeed, we ran some simulations with 
2 and 3 filaments and with different accretion rates for 
each filament: there is still a strong lopsidedness, that can 
reach values of < A\ > larger than 0.10. We show for in¬ 
stance one case, that fits the morphology of NGC 1637, in 

Fig. El 


When accretion is stopped in the simulation, the lop¬ 
sidedness has a typical lifetime of 3 Gyrs (see Fig. ITTIl . It is 
thus rather long-lived, but its lifetime seems a bit smaller 
than in the minor merger and distant interaction scenar¬ 
ios. This is maybe due to the smaller perturbations of the 
halo in the case of gas accretion, while the halo can help to 
make anm = l mode long-lived (e.g., Jog 1997). However, 
a detailed understanding of the lifetime of lopsidedness in 
each scenario is beyond the scope of this paper. 

Lopsidedness and other properties Gas accretion can be 
the clue to explaining strong lopsidedness in galaxies that 
have no sign of recent interaction/merger. We have shown 
that disks do not remain strongly lopsided long enough for 
these situations to be explained by mergers or interactions 
that have occurred long ago. With the accretion of gas, 
we can explain such cases, as for instance NGC 1637, the 
morphology of which is well reproduced in one of our sim¬ 
ulations with accretion along 2 filaments (see Fig. I I81) . If 
asymmetrical accretion of gas with high rates is frequent, 
strong lopsidedness in isolated galaxies can be frequent, 
too. 

The relation between A\ and A 2 is also a natural con¬ 
sequence of gas accretion, for this favors both lopsided¬ 
ness and m = 2 arms/bars at the same time (Bournaud & 
Combes 2002). This scenario can thus explain why disks 
with strong m = 2 asymmetries are more lopsided. 

Gas accretion can also explain why there is a higher 
fraction of lopsidedness and stronger lopsidedness in late- 
type galaxies. First, late-type galaxies are generally less 
massive than early-type, so if they accrete gas at the same 
absolute rate, this can represent a large relative amount of 
mass and trigger a larger lopsidedness. Second, late-type 
galaxies contain more gas, which could be due to the fact 
that they have accreted gas more recently (but there are 
other possible explanations). Finally, gas accretion tends 
to make a galaxy evolve towards late-type (Bournaud & 
Combes 2002), so the most lopsided galaxies in this sce¬ 
nario would naturally tend to be later-type, in contrast to 
what occurs with galaxy mergers that can trigger strong 
lopsidedness but make the galaxy become earlier-type at 
the same time. The correlation between < A\ > and the 
Hubble type T would thus be naturally explained in the 
frame of cosmological gas accretion. 

Then, the observational results that were not explained 
correctly by the interaction and merger scenarios can be 
accounted for by gas accretion. Interactions and mergers 
can still produce lopsidedness in a large number of galax¬ 
ies, but it is required that a large fraction of lopsided disks 
are the result of accretion to explain observational prop¬ 
erties of galaxies. 

This is fully compatible with the high fraction of lop¬ 
sided disks in HI (Sect. 12.511 . Indeed, this observational 
result does not rule out that some lopsidedness can be 
caused by galaxy interactions or mergers, but at least 
it agrees with the accretion scenario, according to which 
gaseous disks should be at least as lopsided as stellar disks. 
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Fig. 18. Snapshot of a simulation of asymmetrical gas accretion (left), giving an overall aspect rather similar to 
NGC 1637 (right: NIR map from the OSUBGS data, after deprojection). A strong lopsidedness is present in this very 
isolated galaxy. We have used two ’’filaments” to provide a more realistic case: the main filament shown in Fig. ED 
is present, with an accretion rate of 4 Mg yr , and a second filament, symmetric of this first one with respect to 
the galaxy center, has an accretion rate of 2 Mg yr -1 . These accretion rates are for a galactic mass of 10 11 Mg, and 
correspond to doubling the mass of the galaxy in about one Hubble time. 


5. Discussion and conclusion 

We measured the Fourier amplitudes of the m = 1 compo¬ 
nent indicating lopsidedness in the stellar disk for a sample 
of 149 galaxies observed in the NIR, 149 of them suitable 
for the study of lopsidedness (inclined by less than 70 de¬ 
grees) . At least two thirds of these galaxies are found to be 
significantly lopsided, in the sense that the m = 1 asym¬ 
metry is too large to simply result from internal mech¬ 
anisms (< Ai > larger than 0.05). The mean value of 
< Ai > in our sample of galaxies is 0.11, and 34% of 
the galaxies show values higher than this. We have shown 
that disks are more lopsided in late-type galaxies, s well 
as when the m = 2 asymmetries (bars and spiral arms) 
are strong. However, the presence of a strongly lopsided 
disk is uncorrelated to the presence of interacting compan¬ 
ions. Many isolated galaxies with no evidence of a recent 
interaction or merger are strongly lopsided. 

Our numerical simulations seek to understand the 
physical origin of lopsidedness in spiral galaxies. We first 
studied distant tidal interactions between galaxies, and 
minor galaxy mergers (not major mergers, for they result 
in elliptical-like systems). The second scenario is much 
more likely. In both situations, strong lopsidedness can 
arise. Yet, both scenarios fail to explain the large fre¬ 
quency of strong lopsidedness in galaxies that are very 
isolated with no companion and no sign of a recent merger. 
They also contradict the observed correlations of the disk 
lopsidedness with the Hubble type, and with the pres¬ 
ence of m = 2 asymmetries. We then simulated an ar¬ 
bitrary model of accretion along one single filament, but 
also more realistic cases with several filaments, which can 
explain both the statistical properties of the sample and 
the morphology of some well-known cases. 

Our study is based not only on the quantitative val¬ 
ues of < Ai >, but also on qualitative relations between 


< Ai > and other parameters. This is important, for the 
quantitative comparison of < Ai > values alone could 
be questioned, as there might be some systematic differ¬ 
ences between the density in simulations and the observed 
’’density” computed from NIR images, especially if the 
mass-to-light ratio of stars varies with radius. However, 
this does not affect the dependence of < Ai > on other 
physical quantities. Both results (values of < Ai > and 
dependence on other parameters) led us to the same con¬ 
clusion, which is therefore robust. 

As for the possibility that disk lopsidedness could sim¬ 
ply result from internal mechanisms, suggested for in¬ 
stance by Masset & Tagger (1997), we do not rule it out; 
but in our simulations, as well as in many other numeri¬ 
cal works in the literature, we have never observed strong 
lopsidedness arising spontaneously in disks. Internal mech¬ 
anisms are found to lead only to values of < Ai > typ¬ 
ically smaller than 0.05, so they are not the clue to the 
origin of the observed disk lopsidedness 2 , since < Ai > is 
larger than 0.05 in most spiral galaxies. Internal mecha¬ 
nisms could only explain the weakest lopsidedness in our 
sample. 

To explain the observations, we show that a large frac¬ 
tion of lopsided disks result from accretion of gas, with 
a typical accretion rate that corresponds to doubling the 
disk mass in about 10 Gyrs, in agreement with the predic¬ 
tions by Block et al. (2002). Observations of lopsidedness 
in HI disks, as traced by the lopsidedness of the global 
HI profiles, agree with this accretion scenario, but NIR 


2 Internal mechanisms are not found to lead to high values 
of < Ai > between 1.5 and 2.5 disk scale-length, but they can 
trigger stronger lopsidedness at smaller radii. They may be 
responsible for strong inner lopsidedness, but cannot account 
for the large-scale lopsidedness studied in this paper 
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observations of stellar disks enable a better quantitative 
study than HI velocity profiles without spatial resolution. 

Large accretion rates of a few MQ-yr -1 along cosmo¬ 
logical filaments may explain both the ubiquity of lopsid¬ 
edness (this paper) and of spiral arms and bars (Bournaud 
& Combes 2002 and Block et al. 2002). In both cases of 
m = 2 and m = 1 asymmetries star formation makes the 
effects of gas accretion visible in the stellar component. 
The results presented here for lopsidedness can even be 
considered as more robust than our previous study based 
on bars. First, deprojection affects bars more than lopsid¬ 
edness. Galactic disks are never perfectly round, even in 
their outer parts, so there is an uncertainty on the disk 
inclination, and the associated error in the deprojection 
will bias the m = 2 asymmetry, but not the m = 1 one; so 
this observational concern is absent in the present study. 
Second, the present conclusions are supported not only by 
the quantitative values of the asymmetry, but also by cor¬ 
relations with other physical parameters, providing several 
independent proofs of gas accretion. Third, lopsidedness 
mainly concerns large-scale outer disks, while bars are cou¬ 
pled with inner gas inflows and nuclear bars. The relative 
resolution of simulations is higher in the outer disk, so 
results on disk lopsidedness are likely to be more robust 
than results concerning bars. Then, we find a robust con¬ 
firmation of the result inferred from the study of bars in 
Block et al. (2002): without denying the role of interac¬ 
tions and mergers in the evolution of galaxies, we here find 
a new evidence that large-scale gas accretion must play a 
major role in the morphological properties of galaxies. 

Appendix A: m = 1 and m = 2 asymmetries in the 
OSUBGS sample 

This appendix gives the values of < A m > for m = 1 and 
2, as defined in Sect. roi for each galaxy in our sample, 
and the inclination angle i used to compute these values. 
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Table A.l. Inclination and values of < Ai > and < A 2 > 
for galaxies in our sample. 
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